Low-energy excitations of a boson pair in a double-well trap 
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I. INTRODUCTION 

Ensembles of ultracold, trapped atoms provide an ideal 
test system for the study of fundamental quantum prin- 
ciples. The manipulation of atoms with photons, pQ, 
has given rise to the experimental realization of Bose- 
Einstein condensation (BEC) [2 , 3, 4 and, more recently, 
the trapping and manipulation of condensates using op- 
tical lattice potentials [51 [6] [7J [8]. The weak coupling 
of neutral atoms to their environment mean that this 
system of cold neutral atoms, confined by a periodic po- 
tential, may prove useful in the investigation of primitive 
quantum information processing [9]. Indeed, such sys- 
tems have already been used to carry out a two-qubit en- 
tangling operation p~Qj[TT], thereby realizing the crucial 
CNOT gate. At the same time, the spatially periodic na- 
ture of the system makes it ideal for the detailed study of 
solid-state Hamiltonians [T2j [13j [14] . The benefit of this 
artificial system, in this regard, lies in the fact that the 
experimentalist can easily vary external control parame- 
ters (e.g. laser intensity or wavelength), thereby varying 
particular parameters of the system Hamiltonian. A de- 
gree of control that is not generally afforded to typical 
solid-state systems. 

The dynamics of a system of ultracold atoms, con- 
fined by an optical lattice potential, can be accurately de- 
scribed within the framework of the Bose-Hubbard model 
[12j [13] . In this model the system Hamiltonian is param- 
eterized by the tunnelling strength between adjacent lat- 
tice sites, J, and the on-site interaction energy, U. The 
Hamiltonian describing the system dynamics can then be 
written as 

H = Jj^blbj + ^>n, + t/^n^-l) , (1) 

where is the annihilation (creation) operator for an 
atom at the lattice site i and hi = b\bi is the number 
operator for that site. Parameter q is the single-particle 
energy at lattice site i and will vary with i for an inhomo- 
geneous lattice. Implicit in this model is the assumption 



* Electronic address: j.f.mccann@qub.ac.uk 



that the dynamics of the system is dominated by single- 
and two-particle effects. In this way, the system of two, 
confined, interacting particles represents the fundamen- 
tal building block for the understanding of these many- 
body systems. Furthermore, continual advancement in 
optical lattice technology means that it has become pos- 
sible to confine small numbers of atoms (e.g. 1 or 2) at 
individual lattice sites, effectively realizing a system of 
two trapped atoms. 

For low-energy collisions the particle interactions can 
be accurately represented within the pseudopotential ap- 
proximation [15 . The eigenstates for a system of two 
particles, interacting via a pseudopotential, can be deter- 
mined analytically for both isotropic [16] and anisotropic 
[T7] harmonic traps. Under such confinement, the 'free- 
space' pseudopotential approximation is found to be suf- 
ficiently accurate provided the length scale associated 
with the particle-particle interactions (a) is short com- 
pared to the length scale of the confining potential (L) 
[T8] , For the case in which a and L are compara- 
ble, one may introduce an energy-dependent scattering 
length and solve for the eigenenergies of the system self- 
consistently [H [20] [21] . 

In addition to providing small numbers of particles at 
individual lattice sites, optical lattice experiments also 
allow for the realization of quasi-one and -two dimen- 
sional systems [22, 23 . Simply increasing the confining 
potential steeply in one or two of the transverse directions 
will effectively 'freeze out' the corresponding degrees of 
freedom [24j [25] . Such systems of reduced dimensional- 
ity can also be achieved using optical or magnetic atom 
waveguides. The theoretical treatment of the particle- 
particle interactions in such low-dimensional geometries 
has been previously considered. For a quasi-one dimen- 
sional (quasi- ID) system it was found that the scatter- 
ing could be treated in terms of a ID, zero-ranged 5- 
potential, renormalized according to the confining poten- 
tial [26]. The physical realization of such quasi- ID trap 
geometries and recent advances in the tuning of atomic 
interactions using Feshbach resonances have permitted 
the study of previously inaccessible regimes. Notably, the 
ID system of impenetrable bosons, or so-called Tonks- 
Girardeau gas [27, 28 , has commanded considerable ex- 
perimental [22j[23] and theoretical [29} l3Ql I3T1 l32l l33l 134] 
interest in recent years. 
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In [34] the detailed theoretical study of two interacting 
particles in a (5-split harmonic potential was considered. 
The DVR techniques, [35] [36], employed in [34 to study 
the £-split trap potential can be easily adapted to other 
types of confining potentials. In the current article we 
utilize these same numerical techniques to study a proto- 
typical two- well trap, defined by V (x) = A [x 4 — kx 2 ] . 
The eigenspectrum for this two-particle system is studied 
and properties of the ground and lowest excited states are 
investigated for varying of the barrier height (dictated by 
n) and the strength of the particle-particle interactions. 
Particular consideration is given to the similarities ob- 
served between the the ground state structure in this 
prototypical two-well potential and that of the 5-split 
potential [34] . 

Similar numerical studies of ultracold few-boson sys- 
tems have been recently reported [37] [38j [39]. In 
this work the authors use a multi-configurational time- 
dependent Hartree (MCTDH) method to study systems 
of several bosons in a double- well trap, with narrow width 
Gaussians used to represent both the central splitting 
potential and the interparticle potential. Where com- 
parison is possible, the results of this numerical MCTDH 
study demonstrate qualitative similarity to the results of 
the present study. 

The remainder of this paper is organized as follows. 
In Sec. [Illwe outline the Hamiltonian that shall be con- 
sidered, for two particles confined by a quasi- ID double 
well potential. In Sec. Ill we present the energy level 



spectrum for the single and two-particle systems, illus- 
trating how the spectrum is influenced by barrier height 
and interaction strength. In Sec. IV we examine various 



properties of the two-particle ground state. The prop- 
erties considered include the g round state wavefunction, 
momentum d istribu tions (Sec. IV A[ ) and von Neumann 
entropy (Sec. IV B). Particular emphasis is given to how 
these properties may be influenced by varying the 'exper- 
imentally controllable' parameters of barrier height and 
interaction strength. In Sec. [V]we systematically exam- 
ine these same properties for the lowest excited states 
of this system. Finally, in Sec. [VI 



we summarize our 



findings and make some concluding remarks. 



II. SYSTEM HAMILTONIAN 

Consider a system of two interacting particles confined 
in two dimensions by means of a 'tight' harmonic poten- 
tial, having trapping frequency co_\_ and associated length 
scale d± = ^/h/mcj±. In the remaining third dimension, 
the confining potential is, relatively, 'loose' and has the 
form 



V (x) = A [a 



KX 



(2) 



The parameters A and k determine the precise form of 
the double-well potential. It is straightforward to verify 
that the two minima of this double-well potential are 
located at x m in = ±^ac/2, with the potential at these 



minima being V (x m in) = —An 2 /A. The well separation, 
x m i n , and the barrier height, V(x m in), are controlled by 
the parameter n. 

As a result of the large energy level separation, asso- 
ciated with the transverse eigenstates (fuj±), the trans- 
verse motion of the particles is 'frozen out'. In this way 
the particles are confined to the lowest motional state 
in each of these transverse directions. In this case the 
system is quasi- ID and may be effectively described by 



H 



E 

i=l,2 



"2m dx} + (Xi 



+g 1D 5 (x 2 - xx) 



(3) 

Here, m is the mass, and x\ and X2 are the coordinates of 
atoms 1 and 2, respectively. The quantity gm represents 
the particle-particle interaction strength, and is related 
to the ID s-wave scattering length (am) through gm = 
—2h 2 /ma\D • In turn, am is related to the 3D s-wave 
scattering length, a^o, through am = —d 2 L /2asD(l — 
CasD/dj_) , where C is a constant and has approximate 
value C= 1.4603 [26]. 

In the limit of tight confinement the free-space pseu- 
dopotential approximation, for the particle-particle inter- 
actions, becomes compromised p~8j[l9]. ^ n ^is case > one 
may obtain the eigenenergies for the system by employ- 
ing an energy-dependent scattering length and solving for 
the energy eigenvalues self-consistently [20j [21] [40] . For 
current purposes it is supposed that we are in the regime 
for which the pseudopotential approximation is still valid 
and the ID collisional coupling, gm, ac ts as a parameter 
for the system. 

The aim is to study the eigenvalues and eigenvectors 
for the 2D Hamiltonian given in Eq. (J3|. To facilitate 
this we introduce the scaling; for i = 1, 2. Under 

this rescaling the time-independent Schrodinger equation 
(TISE) can be written as 



HVi (x u x 2 ) =E^i (x u x 2 ) 



where 



H 



ma* 



= 1,2 



1 d 2 

23x1 



I X a HvX „■ J 



(4) 

+g 1D 5 (x 2 - xi) 
(5) 



Here the scaling factor, a, has been chosen such that 
Ama 6 



h 2 



1 . 



Consequently, 



— / Am \ 1/3 

K vw) K ' 
9iD = w [jL) 9id and 



E 



K 2 



\AmJ 



(6) 



E . 



For convenience we shall drop 'bar' on all quantities and 
use, exclusively, the scaled quantities just described. 
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III. 



ENERGY SPECTRUM 



The eigenspectra for the single- and two-particle sys- 
tem are obtained, subject to the scaling introduced in 
the previous section. A cartesian DVR [35] is used to 
discretize the spatial coordinates x\ and x 2l see [34] for 
details. The discretization scheme used in these calcu- 
lations employs N = 61 mesh points in each dimension 
with a mesh spacing of h — 0.16. Consideration is mainly 
limited to the four lowest eigenvalues of the two-particle 
system (i.e. the lowest band), in particular we are inter- 
ested in the behavior as the parameters k and guj are 
varied. 



A. Single-particle spectrum 

Subject to the scaling introduced in Sec. [ll]the TISE 
for the single-particle system is simply 



1 d 2 

2dx~^ 



+ (x 4 



nx 



_ ^single 



(7) 

The single-particle spectrum is presented in Fig. [I] One 
can see that as the parameter k is increased the lowest 
eigenvalues are pulled downwards in energy as V (x m in) 
becomes increasingly negative. At the same time one 
observes the degeneracy of energy levels as n is increased. 




FIG. 1: Lowest energy eigenvalues for a single particle in 
the potential V (x) — x A — kx 2 . The eigenstates alternate 
between states of odd parity (dashed lines) and even parity 
(solid lines). As the barrier is introduced these states pair up. 



the two-particle system by \E^, the eigenstates for the 
lowest band are then denoted by \I/o, ^1, ^2 and ^3 (see 
Sec. |V| for details). This lowest band corresponds to 
the four lowest levels in Fig. [2|a), representing the two- 
particle system in the absence of interactions. In this 
non-interacting regime the two-particle eigenstates, un- 
der exchange symmetry, are 

^o 1 («; x u x 2 ) =u («; x x ) u («; x 2 ) , 

±Ui (K]Xi)Uo (k.',x 2 )] , 
(ac; xi,x 2 ) =ui (ac; x±) U\ (ac; x 2 ) , (8) 

with the two-particle eigenenergies given by corre- 
sponding combinations of the single-particle energies, 
Ef ngle (K). From Fig. ^a), one notes that as k is in- 
creased all two-particle eigenstates in the lowest band 
become degenerate. This degeneracy follows automat- 
ically from the degeneracy of the states uq (ac; x) and 
u\ (k;x) seen in the single-particle case (see Fig. [I]). 
The states, which are symmetric (solid lines) and anti- 
symmetric (dashed lines) under exchange, are indicated, 
corresponding to boson and fermion pairs. 

The effect of introducing interactions between the two 
bosons is displayed in Fig. |2|b) - (d). In Fig. [2^b) a 
scaled interaction coupling of guj = 1 is considered. The 
symmetric states are shifted upwards in energy as a re- 
sult of the repulsive interactions while the antisymmetric 
states remain unaltered. In the limit of large k one now 
observes two pairs of degenerate levels, as opposed to the 
set of four degenerate states seen in the non-interacting 
case. The energy separation of these two pairs of levels 
is monotonically increasing with increasing k. Increas- 
ing the interaction coupling further, Fig. [2jc) , leads to 
one of the symmetric states being promoted above the 
higher-lying antisymmetric state for small k, but with 
increased k the normal ordering is restored. Finally, Fig. 
|2|d) depicts the same spectrum in the limit of strong 
repulsion: gin = 10. The lowest symmetric state now 
follows closely the energy profile of the lowest antisym- 
metric state. This feature is a universal property for a 
system of strongly interacting bosonic particles in ID. In 
the limit of g\r> — > 00 the bosonic particles become im- 
penetrable, and one enters the so-called Tonks-Girardeau 
regime [27, 28 . The Fermi-Bose mapping [28, 29 allows, 
for example, the ground state of the bosonic system to 
be given by 



B. Two-particle spectrum 

Extending consideration to the two-particle spectrum, 
we focus attention on the two-particle eigenstates belong- 
ing to the lowest band. Denoting the i th eigenstate of 



^ (x u x 2 ) = \V?(x u x 2 )\ . (9) 

The similarity of the energy of the symmetric ground 
state and the lowest antisymmetric state, seen in Fig. 
|2jd), is an indication that the system is approaching this 
Tonks-Girardeau regime. 
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FIG. 2: Lowest energy eigenvalues for the system of particles in a double- well potential of the form V (x) = x A — ftx 2 , as 
a function of the barrier height, ft. The spectrum is displayed for four different values of interaction strength: gm = (a), 
giD = 1 (b), giD = 2 (c) and gm = 10 (d). Introduction of the interaction coupling has the effect of shifting the symmetric 
states (solid lines) upwards in energy, whilst the antisymmetric states (dashed lines) remain unaffected. 



IV. GROUND-STATE PROPERTIES 

The ground state wavefunction, \I/o (#1? #2)? f° r two in- 
teracting particles, is presented in Fig. [3j The individual 
color scale plots will be referenced using standard (row, 
column) matrix notation. 

In the non-interacting case and for ft = 0, plot (1,1), 
the wavefunction is distributed fairly isotropically about 
the centre of the trap. Moving down this column, in- 
creasing ft, the wavefunction expands slightly in both di- 
mensions and takes on a more rectilinear appearance, e.g. 
plot (3,1). For small values of ft, the potential resembles 
that of a square well. The x 4 term gives rise to a steep 
boundary and the distribution of the two (independent) 
particles will be quite uniform, leading to the distribution 
seen in (3,1). As the value of ft is increased the wavefunc- 
tion begins to segregate into four quadrants with sup- 
pression in the region of the barrier (i.e. along the lines 
x\ = and x^ =0). This effect is seen, quite markedly, in 
plot (4,1). Considering the energy level spectrum in Fig. 
|2|a) one can see that for ft = 2 one has not yet reached 
the insulator limit, whereas for ft = 5 one is deep within 
this insulator regime, for which degeneracy is observed 
for the four lowest two-particle levels. 

Turning to the fourth column of Fig. |3j plot (1,4) 
shows the case of no barrier (ft = 0). The repulsive inter- 



action precludes any overlap of the particles. The effect 
of increasing the barrier height to ft = 1, ft = 2 and 
ft = 5 can be seen in plots (4,2), (4,3) and (4,4), respec- 
tively. Again, for small values of ft, the wavefunction 
distribution expands slightly in {x\,x<i) space, but now 
the presence of repulsive interactions distorts the wave- 
function along the line x\ = — x<i- In the insulator limit, 
as we have for (4,4), one sees that the wavefunction has 
split into two clear lobes. 

The behavior observed in Fig. [3] correlates closely to 
the behavior reported for the 5-split potential in [34] : 
the segregation of the wavefunction distribution into four 
quadrants, and the vacancy of two of these quadrants 
owing to the introduction of repulsive interactions. These 
features are essentially generic for double-well systems. 



A. Momentum distribution 

The reduced single-particle density matrix (RSPDM) 
has proven to be an extremely useful mathematical con- 
struct in the analysis of pair correlations [41]. For 
the two-particle system considered here, the RSPDM, 
Pi(x,x r ), for a given eigenstate, \£i (#1, #2), is defined 
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31D = jJIC' = 1 31 D = 2 gin = 5 
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ft = 1 



Kr = 2 




FIG. 3: Color scale plots of the ground state wavef unction, 
(xi,X2), for a boson pair in a double-well potential. The 
color scale runs from blue (largest negative value) through to 
red (largest positive value). The different columns represent 
different values for the interaction coupling. The values con- 
sidered are gm = (column 1), gm = 1 (column 2), gm = 2 
(column 3) and gm — 5 (column 4). In each case the effect of 
varying the barrier height is illustrated down a given column. 
Row 1 corresponds to k = 0, row 2 to k = 1, row 3 to k = 2 
and row 4 to k = 5. Results have been obtained from the 
DVR method with N = 81 mesh points in each dimension 
and a mesh spacing of h = 0.14. Each individual plot spans 
the range —5.6 < x±,X2 < 5.6. 



to be 



Pi 



/+oo 
Vi(x,X 2 )Vi(x',X2)dX2 • (10) 
-oo 



This object has been analyzed in detail for the ground 
state of two particles in a <5-split potential [34 . The be- 
havior of po (x, x') for the double-well, presented here, 
exhibits the same gross features as have been observed 
in [34] for the 5-split trap problem. Instead, in this sec- 
tion we focus on the momentum distributions for this 
system. 

The reciprocal momentum distribution for the i th 
eigenstate, rii (fc), is calculated from the corresponding 
reduced single-particle density, pi(x,x'), through Fourier 
transform 



Pi (x,x')e %k ( x x> )dxdx' 



(11) 

where f_™ ni (k) dk = 1 . Equivalently, one may obtain 
the momentum distribution for this eigenstate by con- 
sidering the diagonalization of pi (x,x'). Specifically, the 



eigenvalue equation is 

Pi (x, x') (j>ij (x r ) dx' = Xijfcj (x) 



f 

j — c 



(12) 



where A^, represents the fractional population of the 
'natural orbital' <pij (x) such that ^2 ■ \j = 1, for each 



i. Using numerical quadrature allows one to rewrite (12) 
as a linear equation. The momentum distribution, rii (fcj, 
may then be obtained from the relation 



(13) 



where pij (k) denotes the Fourier transform of the natural 
orbital <pij (x), 



Pij (k) 



1 



+oo 



<t>ij (x) e 



//>•./ 



dx 



(14) 



The momentum distribution for the ground state is pre- 
sented in Fig. [4] The distributions presented correspond 
to g 1D = (a), g 1D = 1 (b), g 1D 2 (c) and g 1D 5 
(d). Also within each figure, the distributions arising for 
several different values for the barrier height (k) are il- 
lustrated. In the non-interacting case, Fig. [4^a), one ob- 
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FIG. 4: Momentum distribution for the ground state of a 
system of two bosons confined by a double- well potential. 
Four different values of interaction strength are considered: 
(a) giD = 0, (b) g 1D = 1, (c) gm = 2 and gm = 5. In 
each figure the effect of varying the barrier height is also il- 
lustrated. Barrier heights considered are k = (solid line), 
k = 1 (dashed line), k = 2 (dot-dash line) and k = 5 (dotted 
line) . 



serves an initial peaked distribution for n = 0. Increas- 
ing k enhances the peak and narrows the distribution. 
With increasing barrier, the ground-state wavefunction 
adapts to spread over the available interval, leading to 
this reciprocal narrowing in momentum space. Further 
increase in n means that the particles begin to experience 
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the effect of the double-well. As the particles are non- 
interacting, the system displays a single-particle behav- 
ior. For a value of ft = 5 (insulator regime) the particle 
splits between the wells and the momentum distribution 
displays prominent second-order peaks, seen in Fig. |4ja). 
The momentum distributions can be observed by scatter- 
ing or free expansion of the particles in the absence of a 
confining potential. From this perspective the second- 
order peaks correspond to the interference fringes that 
arise from two coherent matter wave sources. 

Introducing an interaction encourages localization and 
has the effect of removing these secondary peaks. Once 
again, for small values of ft, e.g. ft = 1 (dashed line) 
and ft = 2 (dot-dash line), the momentum distribution 
becomes increasingly peaked and narrower. In the pres- 
ence of interactions the particles are restricted to separate 
wells and the interference effects are lost. In addition, the 
localization of the particles leads to a broadening of the 
momentum distribution, as observed in Fig. |4jb), (c) and 
(d), for ft = 5 (dotted line). It is also observed that, in 
the absence of any barrier, ft = (solid line), one sees the 
emergence of higher-energy wings for increasing interac- 
tion, guj. Similar high-energy wings have been reported 
in the TG regime for free space, [26], and harmonic con- 
finement, [3T] . 



B. Von Neumann entropy 

Entanglement is a fundamental expression of informa- 
tion content and is responsible for the increased efficiency 
of some quantum algorithms over their classical counter- 
parts. Previous authors have shown that the von Neu- 
mann entropy of the RSPDM is a good measure of entan- 
glement for a system of two bosons [42 , 43 , 44 . For the 
case of two indistinguishable particles, determination of 
whether or not the two subsystems are entangled requires 
that one considers both the von Neumann entropy of the 
reduced single-particle density matrix, and the Schmidt 
number i.e. number of non-zero eigenvalues (A^) ob- 
tained in the diagonalization of Eq. (12), [45 , 46] [47]. 



Examining the lowest solid line (ft = 0), when no in- 



We use the von Neumann entropy to quantify the entan- 
glement in the position coordinates, x\ and x 2: of the 
particle pair. 

Following the diagonalization of the reduced single- 
particle density, Eq. ( [12] ), the von Neumann entropy 
for the i th eigenstate of the two-particle system (Si) is 
obtained from, 



Si 



13 l°g2 ^ij 



(15) 



1. Variation of von Neumann entropy with interaction 
strength 

Variation of the von Neumann entropy with Q\d, for 
the ground state of this system, is plotted in Fig. [5] 
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FIG. 5: Von Neumann entropy, S, as a function of the inter- 
action strength, gm- This dependence is also illustrated for a 
number of different barrier heights: ft = (solid line), ft = 1 
(dashed line), ft = 2 (dot-dash line), ft = 3 (dot-dot-dash line) 
and ft = 4 (dotted line) . One observes an increased sensitivity 
to gm as the barrier is raised. 



teractions are present (giD — 0) then S = as one 
expects. The product states (with correct symmetriza- 
tion) given in Eq. ([8| represent the eigenstates of the 
non-interacting system. Introducing a small interaction 
has the effect of introducing correlations and results in 
a non-zero entropy. Increasing the interaction strength 
leads to an increasing entropy, saturating at S « 1, 
as for the harmonic potential [34] [44]. This behavior 
can be related to fermionization. As the repulsive in- 
teractions increase, the system enters the TG regime. 
In this regime the ground state of the system can be 
represented by the corresponding system of two non- 
interacting fermions, with correct symmetrization. In 
terms of the eigenfunctions prescribed in Eq. (J8|, the 
ground state of the system is given by l^/™ (k]Xi,x 2 )\ = 

^= [u (ft; xi) ui (ft; x 2 ) - ui (ft; x{) u (ft; x 2 )] • The an- 
tisymmetric state, ^1, in the presence of point-like in- 
teractions, will always give 5 = 1. The ground state of 
the system becomes degenerate with this antisymmetric 
state in the limit of hard-core interactions. 

Fig. [5] also displays the effect of increasing the bar- 
rier height, ft. As the system tends towards the insulator 
regime, the entropy of the system becomes increasingly 
sensitive to changes in g\D, about g\r> = 0. This effect 
was also reported in [34] for the 5-split trap, suggest- 
ing that this is another generic feature associated with 
double- well potentials. The increased barrier height re- 
duces tunneling between the wells. For any increase in 
the interaction strength, the two-particle wavefunction 
will attempt to redistribute so as to minimize this inter- 
action. However, with the increased barrier height the 
wavefunction is forced to remain more localized, and is 
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restricted in its redistribution. 



2. Variation of von Neumann entropy with barrier height 



separation of wells only becomes apparent for large values 
of ft (i.e. k > 3). From Fig. [6] one can identify 2 < k < 3 
as the interval over which the entropy makes its most 
rapid variation. 



Variation of the von Neumann entropy with ft, for the 
ground state of this system, is plotted in Fig. [6j The 
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FIG. 6: Von Neumann entropy, S, as a function of the bar- 
rier height, ft. The dependence is examined for a number of 
different interaction coupling strengths: gm = 1 (solid line), 
gm = 2 (dashed line), gm = 5 (dot-dash line) and gm — 10 
(dotted line). In all cases S saturates at a value of unity in 
the limit ft — > oo. 



basic trends bear a striking resemblance to those seen for 
the (5-split trap, [34]. Specifically, one observes that the 
initial entropy of the system (i.e. for ft = 0) is dictated by 
the interaction strength of the system, g 1 ^. The larger 
is gin, the larger is the initial value of 5, as is consistent 
with Fig. [5] Increasing the height of the barrier then has 
the effect of increasing the entropy of the system towards 
5 = 1. In the limit of large barrier heights the entropy 
of the system saturates at S = 1, regardless of the value 
of interaction strength (the notable exception being the 
non-interacting case, for which S is identically equal to 
zero for all ft). This saturation at S = 1 corresponds to 
the loss of entanglement. 

When the system enters the insulator limit there is an 
implicit exchange uncertainty in the state of the system, 
arising from the indistinguishable nature of the particles. 
As such, these correlations cannot be exploited in any 
meaningful quantum information protocol and the sys- 
tem is regarded as non-entangled. This diagnosis also 
follows from the criteria set out in [45j [46] [47] as, in the 
limit ft — > 00 then 5 — > 1 and the Schmidt number can 
be seen to approach a value of 2 (not shown here). By 
the criteria outlined in [46], any state for which S = 1 
and with a Schmidt number of 2 must be regarded as 
non-entangled. 

In contrast to the <5-split trap, the entropy dependence 
S (ft) for the double-well system is quite sigmoidal. The 



3. Von Neumann entropy in the Bose- Hubbard model 

One may examine the von Neumann entropy of the 
ground state within the formalism of the Bose-Hubbard 
model presented in Eq. ([!]). Using a Fock basis for the 
two-particle system of the form Itiltir), where n L ^ rep- 
resents the number of particles in the left (right) well, 
leads to three basis states : |20), |11) and |02). 

Thus, in terms of this basis the Hamiltonian ([!]) may 
be written in matrix form as 



H 



/2e + 2U V2J \ 

V2J 2e V2J 
\ V2J 2e + 2U J 



(16) 



with eigenvalues of this Hamiltonian follow from some 
simple algebra 



E_ =2e + U- \/U 2 + 4 J 2 
E mid =2e + 2U 



E+ =2e + U + \/U 2 + 4 J 2 . 



(17) 



The eigenvectors corresponding to the eigenvalues pre- 
sented in Eq. (17) are found to be 



mid 



) = 



V2 



|*±>=A/- ± 



1 

2V2J 



U±VU 2 +4J 2 

1 



(18) 

where A/± represent normalization factors and l^mid) 
has odd inversion symmetry. 

Given the two-particle ground state in the Fock basis, 
, one may determine the reduced single-particle 
density matrix by tracing over the degrees of freedom of 
either particle. The single-particle basis states may be 
represented in the form \til,tir) as 1 10 ) and 1 01 ) . In 
turn, the two-particle basis states can be written in the 
symmetric form: 



|20> = |10> 1 ® |10> 2 

|11)=-^[|10) 1 ® |01) 2 + ioik® |10) 2 ] 



|02) = 101)! ® |01) 2 



(19) 
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The eigenvalues of this RSPDM are found to be 
4J 2 4J 



Ai =A/ r2 . 
A 2 =Ml 



1 



(U-y/U 2 + 4J*) c/-VC/ 2 + 4J2 

4J 2 4J 
([/ _ + 4 + J7 - ^^+4/2" 



(20) 



The variation of the ground-state entropy (S) with the 
model parameters J and U is depicted as a surface plot 
in Fig. [7| It is noted that the qualitative behavior of the 




FIG. 7: Von Neumann entropy of the ground state as a 
function of the Hamiltonian parameters J and U. The limit 
J/U — > corresponds to the insulator limit and U/J — > to 
the conductor limit. 



entropy displayed in Figs. [5] and [6] for varying girj and ft, 
respectively, is reflected in the Bose-Hubbard model with 
variation of the parameters U and J. This tight-binding 
approximation is poor in the limit U/J — > 0, but is an 
accurate representation in the insulator limit. 



V. EXCITED STATES 

Attention is now turned to the three lowest excited 
states which, together with the ground state, represent 
the lowest energy band of the two-particle double-well 
system. In the non-interacting case, with spectrum de- 
picted in Fig. [2|a), the three lowest excited states may be 
represented as given in Eq. (|8|, with one of these states 
being antisymmetric and two of them symmetric. Varia- 
tion of parameters n and gm can lead to reordering of the 
energy eigenvalues, as observed in Fig. [2] However, in 
this section, the study of the excited states of this two- 
particle system will be restricted to these three states 
of the lowest band, identifiable through their symme- 
try. Henceforth the term 'first-excited state' refers to the 
lowest energy antisymmetric state (^i), 'second-excited 
state' refers to the second-lowest energy symmetric state 



(^2) and 'third-excited state' refers to the third- lowest 
lying symmetric state (^3). 



A. Two-particle excitations 

The wavefunctions for ^1,2,3 are represented, by means 
of color scale plots, in Figs. [5J [9] and [TUJ respectively. 
Again, the standard (row, column) notation is used to 
reference individual subplots. The color scale is consis- 
tent across all wavefunction plots, permitting direct com- 



parison between Figs, pjpl^land 10 
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FIG. 8: Color scale plots of the first-excited eigenfunction, 
\£i (xi,X2), for a pair of particles in a double-well potential. 
The color scale runs from blue (largest negative value) to red 
(largest positive value). The different columns represent dif- 
ferent values for the interaction coupling. The values consid- 
ered are gm — (column 1), gm — 1 (column 2), gm — 2 
(column 3) and gm = 5 (column 4). In each case the effect 
of varying the barrier height is illustrated down a given col- 
umn. Row 1 corresponds to k = 0, row 2 to k = 1, row 3 
to k = 2 and row 4 to k = 5. Notice how the antisymmetric 
state is completely independent of the interaction parameter, 
gm- Results have been obtained using the same DVR method 
described for Fig. [3] 



Fig. [8] represents the ground state for a system of 
two spin-aligned fermions (^1), which is identically zero 
along the line x\ = £2 and, thereby, unaffected by the 
zero-ranged interaction. Considering Fig. [§J moving 
along a given row (i.e. increasing repulsion for a fixed 
barrier), the wavefunction plots remain unchanged, illus- 
trating the independence of this state with respect to 
interaction strength, gi£>. As k increases (i.e. down 
any column) the positive and negative lobes along the 
x\ = — X2 diagonal become more widely separated indi- 
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eating isolation into separate wells. Once again the wave- 
function density in these two quadrants correspond to the 
situation where particle 1 is in the left well (x\ < 0) and 
particle 2 is in the right (x2 > 0), and vice versa. In the 
limit of a large barrier, the ground state becomes degen- 
erate with this antisymmetric state. The wavefunction 
plots are almost identical (except for sign) for n = 5 as 
seen, for example, by comparing Fig. [| (4,4) and Fig. 
[8] (4,4). Furthermore, as already discussed, one expects 
the ground state to become degenerate with this anti- 
symmetric state in the limit of guj — > oo, for all n. This 
degeneracy is evidenced by comparing the fourth column 
in Fig. [3] to any column in Fig. [8] Even at this finite 
interaction strength (gip = 5) the equivalence of these 
two states is apparent. Finally, from each of the plots 
in Fig. [8] it is clear that this eigenstate is of odd parity, 
such that ^1 (#i, X2) = (— #1, — X2). 



gin = 
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FIG. 9: Color scale plots of the second-excited eigenfimction, 
^2 (xi,X2), for a pair of particles in a double-well potential. 
The color scale runs from blue (largest negative value) to red 
(largest positive value). The different columns represent dif- 
ferent values for the interaction coupling and the effect of 
varying the barrier height is illustrated down a given column. 
The same values for the parameters g\r> and k, , as considered 
in Fig. [8] are examined here. Results have been obtained 
using the same DVR method described for Fig. [3] 



Fig. [9] depicts the second-excited state for the sys- 
tem of two bosons in a double-well potential and, as 
with ^1, this state exhibits odd parity: ^2(^1,^2) = 
—^2 ( — #1? ~ x< i)- The case of no barrier (k = 0) and no 
interaction (gm = 0) is illustrated in Fig. [9] (1,1). The 
eigenstate is composed of two lobes which correspond to 
both particles co-existing on the same side of the well. 
In the case of no interactions (giD — 0), illustrated in 
column 1, this symmetric eigenstate is degenerate with 



the antisymmetric state considered in Fig. [8j Repul- 
sive interactions will tend to exclude the wavefunction 
from the line x\ = x^ (e.g. compare (1,1) to (1,3) or 
(1,4)). In the Tonks limit this splits each of the upper 
right and lower left lobes. Considering the effect of the 
barrier in column 3 (gm — 2), the initial wavefunction 
demonstrates the double-lobe structure. As the barrier 
is increased to k = 1, (2,3), and then n = 2, (3,3), the 
wavefunction spreads out in (^1,^2) space. Further in- 
crease of the barrier height causes the system to move 
into the insulator limit, plot (4,3), forming two isolated 
lobes in the upper-right and lower-left quadrants. The 
eigenstate, in this case, corresponds to the physical situ- 
ation of both particles residing in either the left well or 
the right well. 
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FIG. 10: Color scale plots of the third-excited eigenfunction, 
^3 (xi,X2), for a pair of particles in a double- well potential. 
The color scale runs from blue (largest negative value) to red 
(largest positive value). The different columns represent dif- 
ferent values for the interaction coupling and the effect of 
varying the barrier height is illustrated down a given column. 
The same values for the parameters giD and as considered 
in Fig. [SJ are examined here. Results have been obtained 
using; the same DVR method described for Fig. [3] 



Finally, the third-excited state is illustrated in Fig. 
[TO) In contrast to ^1 and ^2, this eigenstate is of even 
parity such that ^3 (£i,£2) = ^3 ( — #1? — x 2j- Scanning 
down column 1: as the system moves into the insulator 
limit, the eigenstate is composed of four equally- weighted 
lobes in the four quadrants, equivalent to the correspond- 
ing ground- state eigenfunction, seen in Fig. [3] (4,1). In 
fact from Fig. [2|a), in the non-interacting case, the four 
lowest eigenstates all become degenerate in the insulator 
limit (k —> 00). As a consequence, the eigenfunction plot 
(4,1) in Figs. |3j |8j [9] and 10 relate to four degenerate 
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states. 

This symmetric eigenstate is non-zero along the line 
x\ — X2- As one increases giD one again observes the 
exclusion of the wavefunction from this line (e.g. exam- 
ining row 1 in Fig. 10). As barrier height is increased the 



wavefunction expands in (#i, #2) space and there is some 
suppression of the wavefunction in the region of the rising 
barrier (i.e. X\ = and x^ =0). In the insulator limit, 
e.g. in plot (4,4) for which k = 5, the wavefunction in 
the off-diagonal quadrants vanishes and one observes two 
double- lobes in the lower- left and upper-right quadrants, 
representing the physical situation where both particles 
reside in the same well. The degeneracy of ^2 and ^3, 
in the limit k — » 00 (seen in Fig. [2| is manifested in the 
corresponding wavefunction plots. This is demonstrated 
by comparing corresponding plots in the bottom rows 
(k = 5) of Figs. [9]and[l0| 

In the limit of large k (and for any positive interaction), 
the ground and first-excited states correspond to the two 
particles in separate wells. By contrast, the second- and 
third-excited states, in the same limit, correspond to two 
particles in the same well. It follows that an increase in 
the repulsive interaction coupling will cause this second 
pair of levels to be shifted upwards in energy. In this way, 
in the k — > 00 limit, one observes the separation of these 
two pairs of levels to increase as g\r> is increased (see 
Fig. [2]). The increasing separation of these levels with 
increasing n is corollary to this. As n is increased the 
particles become more tightly confined to the individual 
wells. This increased confinement, for the upper pair of 
levels, will give rise to an increased interaction of the two 
particles and a subsequent increase in the energy of these 
eigenstates, relative to the lower pair. 



B. Momentum distribution 

The momentum distributions for the excited states are 
calculated as outlined in Sec. IIV Al The calculated dis- 
tributions for the second-excited state (^2? where the 
superscript 'o' indicates the 'odd' inversion symmetry of 
this eigenstate) are displayed in Fig. [II] For k = one 
observes a double-humped distribution that becomes nar- 
rower with increasing n and, in the insulator limit, gives 
way to a single-peak distribution with high-energy tails. 
This is similar to the result for (not shown). An in- 
crease in the interaction coupling has the effect of narrow- 
ing the momentum distribution. Fig. [9] illustrates that 
increasing gm will expand the wavefunction in (£i,£2) 
space, leading to this reciprocal narrowing in momentum 
space. At the same time the increased interaction leads 
to an accentuation of the double-peaked structure, ob- 
served for small n. 

Fig. 12 illustrates the momentum distribution for 
the even-parity state For the non- interacting case, 
9id — (a), a double-mode distribution arises with 
a node at k = 0. This node is accounted for due to 
the separable nature of ^| in the non- interacting limit: 



0.3- 
g0.2 
0.1- 




0.3 

g0.2 
c 

0.1 




1 I 1 I 1 I 1 I 1 I 1 

- ( a ) ^ 9id = ° 


1 | 1 | 1 | ■ | 1 | 1 

" ( b ) S\ 9id = 1 - 

J. Jr*Vl 1 1 1 1 \*L..J 


" (°) /^>- 9lD = 2 

l-.L^ I 1 I 1 1 \.*J... 


- ( d ) /;.9id = 5 - 



-4 -2 2 4 -4 -2 2 4 

k k 

FIG. 11: Momentum distribution for the second-excited state 
(^2) f° r a system of two atoms confined by the double- well 
potential V (x) — (x 4 — kx 2 ). Four different values of inter- 
action strength are considered: (a) gm = 0, (b) gm = 1, (c) 
gm = 2 and gm = 5. In each figure the effect of varying the 
barrier height is also illustrated. Barrier heights considered 
are k = (solid line), k = 1 (dashed line), k = 2 (dot-dash 
line) and k = 5 (dotted line). The momentum distribution is 
identical to that of in the non-interacting limit. Increasing 
the interaction strength (gm) leads to an increasingly peaked 
distribution. 
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FIG. 12: Momentum distribution for the third-excited state 
(^3) for a system of two atoms confined by the double- well 
potential V (x) — (x 4 — kx 2 ). Four different values of inter- 
action strength are considered: (a) gm — 0, (b) gm = 1, (c) 
gm = 2 and gm = 5. In each figure the effect of varying the 
barrier height is also illustrated. Barrier heights considered 
are k = (solid line), k = 1 (dashed line), k — 2 (dot-dash 
line) and k — 5 (dotted line). The separable nature of ^| in 
the non-interacting limit, as seen in Eq. ([8]), leads to node 
in momentum distribution for k — 0. Introduction of inter- 
actions means that ^| is no longer separable and node is no 
longer enforced. 



^cj 1 (ft; £1, £2) = u\ (n\ x\) u\ (k; X2) . Considerable nar- 
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rowing of this distribution is noted as ft is increased 
and in the insulator limit a second pair of smaller peaks 
emerges. This second pair of peaks may be viewed as 
interference fringes from each particle being distributed 
between the two wells - compare Fig. 0(4,1) and Fig. p 
(4,1). 

For increased interaction strength (qid) one contin- 
ues to observe the narrowing of the distribution with 
increased barrier height. However, the presence of the 
interactions causes the node at k = to be removed, as 
one can no longer write the eigenfunction in the separable 
form given in Eq. (J8|. Instead one just observes a strong 
depression of the distribution about k = 0. At the same 
time, the introduction of the interactions has the effect of 
completely removing the double-peaked structure in the 
insulator limit, as is observed for the dotted line (ft = 5) 
in each of Fig. 12 'b), (c) and (d). As seen in Fig. [To! 



in the presence of a finite interaction the wavefunction in 
the off-diagonal quadrants vanishes in the insulator limit, 
and this eigenstate describes a situation where both par- 
ticles occupy one side of the double-well. 



C. Von Neumann entropy 

As for the ground state, one may obtain the von Neu- 
mann entropy for the excited states of the two-particle 
system via diagonalization of the reduced single-particle 
density matrix. In this section the dependence of the 
von Neumann entropy, 5, of the four lowest two-particle 
states, on the interaction strength (gm) and the barrier 
height (ft) is considered. 



1. Variation of von Neumann entropy with interaction 
strength 



Fig. [13] illustrates the dependence of £ on the interac- 
tion coupling, gin (> 0). The dependence is examined 
for four different values of barrier height: ft = (a), 
ft = 2 (b), ft = 4 (c) and ft = 5 (d). For each value of the 
barrier height the entropy of the four lowest eigenstates 
is depicted: ground state (solid line), first-excited state 
(dashed line), second-excited state (dot-dash line) and 
third-excited state (dotted line). The dependence of the 
ground-state entropy on gi^ has already been examined 
in Fig. [5j however it is useful to replicate these plots here 
to help inform the examination of the excited- state plots. 

Several important features are noted. In all cases the 
first-excited state (dashed line) shows no dependence on 
the interaction strength, as is expected owing to the sym- 
metry of this eigenstate. Instead, this eigenstate exhibits 
a value of S = 1 for all guj. This value follows from the 
analytic form for this eigenstate, ^J 1 , given by Eq. ([8]), 
which holds for all values of g\D- At the same time, the 
analytic representations for the three remaining eigen- 
states are also given in Eq. (|8|, for guj — 0. From these 
representations it is clear that, in the non-interacting 



limit, the entropy for the ground state (solid line) and 
third-excited state (dotted line) is always zero, as these 
states may always be represented as direct-product states 
for giD — 0. In a similar way, the second-excited state 
(dot-dash line) always assumes a value of S = 1 in the 
non- interacting limit. Once again, this may be attributed 
to the symmetrized form for this state as given by m 
Eq. (|8]). 



Considering the case of ft = 0, Fig. 13 'a), the ground 
state begins at S = and increases monotonically with 
giD- As giD — > oo one enters the TG regime and this 
ground state (solid line) becomes degenerate with the 
first-excited (dashed line) state and S ~ 1. By contrast, 
the second-excited state begins with S = 1, as discussed, 
and increases with increasing #i£>, but at a much slower 
rate than that exhibited by the ground state. The third 
excited state (dotted line) begins, like the ground state, 
with S = and increases rapidly with increasing inter- 
action strength. 

Increasing the height of the barrier to ft = 2, Fig. 
[l3]fo), one observes qualitatively similar behavior from all 
four states except that each state exhibits a more marked 
variation in S over the range of gw examined. As one 
moves into the insulator regime, e.g. ft = 4, Fig. [l3^c), 
the behavior changes quite significantly. As discussed 
previously, the ground state exhibits a very drastic vari- 
ation with giD, converging very rapidly to S « 1. The 
second-excited state still exhibits the same basic behavior 
as noted for smaller ft but, once again, the increased bar- 
rier height leads to an increased sensitivity of this state 
to variation in g\r>- The third-excited state shows a dis- 
tinct change in behavior for this increased barrier height. 
At small values of interaction coupling (giD < 1) the 
entropy of this state follows closely that of the ground 
state. As interaction strength is increased beyond this 
value then the ground-state entropy begins to plateau 
at S « 1, whilst that of third-excited state continues to 
increase. Increasing the barrier height to ft = 5, Fig. 
[l3^d), moves the system deeper into the insulator limit 
and the behavior demonstrated in (c) becomes even more 
striking. In this case the behavior of the ground-state en- 
tropy is more dramatic, with the entropy saturating at 
S « 1, already, for guj « 0.5. Again the entropy of the 
third-excited state follows this trend identically. How- 
ever, where the entropy of the ground state plateaus at 
5 « 1, the entropy of the third-excited state continues to 
increase and follows now, almost identically, the entropy 
of the second-excited state. A handle on this behavior 
is provided by the wavefunction plots of Figs. [3j [9] and 
[To} One observes that, in this insulator limit (ft = 5), 
the third-excited state, for small g\D, as seen in Fig. [10] 
(4,1), closely resembles the ground state in Fig. [3] (4,1). 
For larger interaction couplings (gin > 1) the eigenfunc- 
tion for this third-excited state, as seen in Fig. [To] (4,2) 
- (4,4), closely resembles that of the second-excited state 
in Fig. |9](4,2) - (4,4). 
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FIG. 13: Von Neumann entropy, 5, as a function of the interaction strength (gm) for the four lowest two-particle states: 
(solid line), (dashed line), ^ (dot-dash line) and \£| (dotted line). The dependence is also illustrated for a number of 
different barrier heights: ft = (a), ft — 2 (b), ft — 4 (c) and ft = 5 (d). 



Variation of von Neumann entropy with barrier height 



The variation of von Neumann entropy with barrier 



height is illustrated in Fig. 14 , for the same four, lowest- 
energy two-particle states. In this case, four different 
values of interaction coupling are presented: guj = 1 
(a), gi D = 2 (b), g 1D = 5 (c) and g 1D = 10 (d). Again, 
in each plot the eigenstates are represented by the same 



line types used in Fig. 13 



Some general features and behaviors can be noted from 
these plots. Again, the first-excited state is observed to 
have an entropy of unity for all gin and k. For k — > oo, 
the entropy of the ground state tends to a value of unity, 
regardless of the value of guj (provided guj > 0). In 
this limit the ground state of the system is described 
by one particle in each half of the double-well potential, 
and corresponds to the Mott-insulator regime. On the 
other hand, the initial value of S (when ft = 0) is sen- 
sitive to gio- The higher the value of g\D, the larger 
is the initial value of S. As S — » 1 in the insulator 
limit, it follows that the entropy of the ground state ex- 
hibits a less dramatic variation with ft, for larger values 
of interaction strength. For all of the symmetric eigen- 
states, i.e. ground (solid line), second-excited (dot-dash 
line) and third-excited (dotted line), as the interaction 
strength is increased the entropy of the eigenstates, in 



the entropy of these symmetric states in the absence of 
a barrier (ft = 0), increases with increasing g\r>. The 
second-excited state (dot-dash line) exhibits an entropy 
that monotonically increases with ft for the range of pa- 
rameter space considered. By contrast, the third-excited 
state (dotted line) exhibits an entropy that both increases 
then decreases with raising of the barrier. 

One will also note that in the limit of large barrier 
heights (i.e. ft — > oo), the entropy of the second- and 
third-excited states tend to the same value. Although 

^d), this fact has also been ver- 
= 10. Once again, a handle on 



not obvious from Fig. 
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general, increases, consistent with Fig. 13 In particular, 



ified for the case of g\r> 
why this happens can be obtained from the wavefunction 
plots for these eigenstates in Figs. [9] and [lOj One can 
see that in the presence of finite interactions, these two 
eigenstates become identical in the insulator limit, except 
for some phase (compare row 4 of these figures). Fig. 14 
also suggests the the value of S to which these two states 
converge, in the insulator limit, is greater than one and 
increases with increasing interaction strength, g\r>. 

This behavior of the entropy may be qualitatively un- 
derstood as follows. Both the second- and third-excited 
states correspond, in the insulator limit, to the phys- 
ical situation of two particles coexisting in either the 
right-well or the left-well. As such, these states may 
be roughly represented by Bel l-type s tates of the form 
l/v / 2(|20)± |02)) - see Sec. |lVB3|for the definition 
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FIG. 14: Von Neumann entropy, S, as a function of the barrier height (k) for the four lowest two-particle states: (solid line), 
(dashed line), ^2 (dot-dash line) and ^| (dotted line). The dependence is examined for a number of different interaction 
coupling strengths: (a) gm = 1, (b) g\D — 2, (c) gm = 5 and (d) gm = 10. 



of these basis states. Such a Bell state carries one e- 
bit of entanglement, with a corresponding von Neumann 
entropy of unity. However, beyond this, there are also 
correlations between the two particles coexisting in the 
same well. As can be seen, for example, from Fig. [10] 
(4,4). Here the repulsive interactions between the parti- 
cles occupying the same well leads to a partition of the 
wavefunction, within each well, into two lobes. Consid- 
ering the double lobe seen in the upper-right quadrant, 
corresponding to both particles co-existing in the right 
well of the double-well. The upper half of the lobe repre- 
sents the situation where particle 1 is on the left of this 
well and particle 2 is on the right, the lower half-lobe cor- 
responds to the reverse of this situation (x\ <-+ X2). In 
this case the correlations in the system are analogous to 
the correlations that are observed for the ground state, 
^o, m the absence of any barrier (k = 0). These cor- 
relations (and therefore S) are seen to increase as the 
interaction coupling is increased. One significant distinc- 
tion exists between these 'single-well' correlations, seen 
in states ^2,3, and the correlations seen in the ground 
state, for k = 0. On increasing the second- and 
third-excited states tend to become more confined and 
the two-particle wavefunction becomes increasingly lo- 
calized in the single-well. However, the particle-particle 
interactions will compete with this effect, attempting to 
keep the two-particle wavefunction spread in space and, 



in particular, minimized along the line x\ = x^. For 
the ground state this particular type of single-well com- 
petition between k and qid is not experienced. So, in 
the insulator limit, the second- and third-excited states 
will have correlations arising from the realization of the 
Bell-type state, and the 'single-well' correlations due to 
the two interacting particles coexisting in the same well. 
This combination of factors leads to an entropy which is 
greater than unity, with the contribution of the 'single- 
well' correlations, in general, increasing with increasing 
interaction strength. 



D. Stimulating two-particle excitations 

The previous results have clearly illustrated that ma- 
nipulations of this two-particle system can be achieved 
through the variation of the control parameters guj and 
k in some adiabatic manner. However, one could also 
consider time-dependent manipulation of the state. Con- 
sidering the insulator limit, one may propose two meth- 
ods of coupling these lowest levels: (a) shaking the trap 
from side-to-side (b) modulating the barrier height (see 
Fig. 15). To first-order, the former represents a dipole 



excitation, capable of coupling |\I/q) and )• The lat- 
ter scheme (to first-order) corresponds to a quadrupole 
excitation, capable of coupling states |\£q) and In 
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this way, by employing such techniques it should prove 
possible to exploit these three lowest eigenstates in order 
to engineer the two-particle state in a time-dependent 
fashion. 
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FIG. 15: Time-dependent manipulation of the two-particle 
state may be achieved by (a) oscillation of the trap from side- 
to-side or (b) the modulation of the barrier height. The for- 
mer process provides a dipole coupling between the ground 
state, l^o)? an d the state 1^2 )■ The second scheme will pro- 
vide a quadrupole coupling between the ground state and the 
state 1^1). (c) shows a schematic energy level representation 
for these processes. 



Further investigation of this idea of time-dependent 
manipulation of the two-particle state could prove a use- 
ful extension to the present study. In particular, a com- 
bination of time-dependent excitation processes and the 
adiabatic variation of control parameters, gin and ft, 
should permit an impressive degree of control over the 
two-particle state, within this system. 



VI. 



SUMMARY 



The system of two interacting particles in a pro- 
totypical double-well potential of the form V (x) = 
A [x 4 — kx 2 ] has been considered. Using a cartesian 
DVR, the eigenspectrum for this system has been studied 
and the four lowest eigenstates have been obtained and 
investigated for varying barrier height and interaction 



strength. For each state the two-particle eigenfunction, 
the momentum distribution and the von Neumann en- 
tropy have been examined. It was found that the ground 
state for this double-well system exhibits behavior that 
closely resembles that observed in a previous study of 
the £-split trap potential, [34 . In particular, the ground- 
state wavefunction is suppressed along the lines x\ = 
and X2 = as barrier height is increased, leading to a 
quadrant separation of the wavefunction. In the presence 
of repulsive interactions (gio > 0) only the contributions 
in the off-diagonal quadrants remain in the insulator limit 
(ft — > 00). In this limit the ground state of the system is 
composed of one particle in each half of the double-well. 
The momentum distributions display an initial narrow- 
ing with increasing barrier height but with a broadening 
and high-energy wings being observed in the insulator 
limit. Furthermore, the secondary peaks observed in the 
momentum distribution for the double well, in the non- 
interacting regime, are quickly suppressed in the pres- 
ence of repulsive interactions. The variation in the von 
Neumann entropy (S) with interaction strength shows 
remarkably similar behavior. In all cases S = in the 
absence of interactions and for gm — > 00, S saturates at 
a value close to unity. Increasing the height of the barrier, 
in each case, has the effect of making the entropy more 
sensitive to changes in the interaction strength, around 
9id — 0. Similarly, the behavior of the entropy with 
varying barrier height exhibits generic features between 
the two double-well systems. In both cases the ground- 
state entropy saturates at a value of unity as ft — > 00, 
regardless of the value of gio- The initial value of S (i.e. 
the value of S for ft = 0) is determined by the strength of 
the interaction, with larger interaction coupling leading 
to larger initial entropy. As such, the sensitivity of S to 
ft is reduced for double-well systems with larger interac- 
tion couplings (<7i£>). This behavior of the ground-state 
entropy is also illustrated within a Bose-Hubbard model, 
wherein the controllable parameters are the on-site inter- 
action (U) and the tunnelling strength (J). 

As well as examining the ground state of this double- 
well system, some of the properties of the three lowest 
excited states have also been studied. Two of these states 
are found to be symmetric whilst one is antisymmetric 
and they constitute the lowest band of the two-particle, 
double-well system. The antisymmetric state is found 
to be completely independent of the interaction parame- 
ter (gin)- However, this state displays a dependence on 
the barrier height and in the limit of a high barrier be- 
comes degenerate with the ground state - corresponding, 
physically, to the situation of each particle residing in a 
separate, isolated well. The von Neumann entropy for 
this antisymmetric state is identically equal to one for all 
ft and gio- 

The second- and third-excited states are symmetric. 
In the insulator limit (provided g\r> > 0) the states be- 
come degenerate and correspond to the physical situa- 
tion where both particles occupy the same well. Both 
eigenstates demonstrate momentum distributions that 
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are double-humped, with the double-hump giving way 
to a single peak in the insulator limit. For gm = 
the third-excited state exhibits secondary peaks in the 
momentum distribution, similar to the ground state in 
the non-interacting regime. The entropy of both states 
increases with #i£>, with that of the third-excited state 
showing a more marked variation. As for the ground 
state, increasing barrier height n has the effect of in- 
creasing the sensitivity of the entropy to variations in 
giD (about gio =0). In the insulator limit the entropy 
of the third excited state is found to follow, almost iden- 
tically, that of the ground state for small gw- As the 
ground state entropy saturates at S « 1, the entropy 
of the third excited state continues to increase and, for 
larger g\D, follows, almost identically, that of the second- 
excited state. Indeed, in the insulator limit and for fixed 
interaction strength, the second- and third-excited states 
are found to have the same entropy (as follows from the 
physical equivalence of these states in this limit). The 
entropy, in this case, is proposed to have two contribu- 
tions due to (i) the realization of a Bell-type state with 
both particles co-occupying either the left or right well 
(ii) single-well correlations, owing to the repulsive inter- 
action of the two particles occupying the same well. 



pared to the £-split trap previously considered. Having 
characterized the properties of the ground and lowest- 
excited eigenstates, the foundation is laid for future in- 
vestigation into state manipulations using this system. 
Future avenues may include the time-dependent manipu- 
lation of states through shaking of the trap, an oscillating 
barrier height or introduction of a constant, or oscillat- 
ing, field gradient. These time-dependent manipulations, 
along with the adiabatic variation of the control param- 
eters gio and ft, should allow for comprehensive state 
engineering within the lowest band of this two-particle 
system. 
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